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Abstract 

We have searched for lepton flavor violating decays of the r lepton with one or two K s 
mesons in the final state. The data used in the search were collected with the CLEO II and 
II. V detectors at the Cornell Electron Storage Ring (CESR) and correspond to an integrated 
luminosity of 13.9 fb^ 1 at the T(4S) resonance. No evidence for signals were found, therefore 
we have set 90% confidence level (C.L.) upper limits on the branching fractions B(t~ — > 
e~K° s ) < 9.1 x 1(T 7 , B(t~ -> n~K%) < 9.5 x 1(T 7 , B(t~ -»■ e~K° s K° s ) < 2.2 x 10~ 6 , and 
B(t~ — ► /j,~K s K s ) < 3.4 x 10 -6 . These represent significantly improved upper limits on 
the two-body decays and first upper limits on the three-body decays. 

PACS numbers: 13.35.Dx, 11.30.Hv, 14.60Fg, 14.40AQ 
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In physics, all fundamental conservation laws are expected to have associated symme- 
tries. Lepton flavor conservation, however, is an experimentally observed phenomena with 
no associated symmetry in the Standard Model. Lepton flavor violation (LFV) is expected 
in many extensions of the Standard Model such as lepto-quark, supersymmetry, superstring, 
left-right symmetric models, and models that include heavy neutral leptons |IJ. Experimen- 
tally, both Super Kamiokande [0] and SNO || observe neutrino oscillation, which may imply 
LFV in the neutrino sector; therefore LFV is expected to occur in charged lepton decay at 
some branching fraction, albeit very small. The r lepton provides a clean laboratory for such 
searches. Ilakovac @] has calculated upper limits on branching fractions for many neutrino- 
less LFV modes within a model involving heavy Dirac neutrinos. The branching fractions 
depend on the heavy neutrino masses and mixings. For the decays t~ — > £~K$ ||, the 
branching fractions are of O(l0~ 16 ), where I can be e or fi. For the decays t~ — > £~KgKg 
or — > £~K + K~, the branching fractions are of O(l0~ 7 ). The decays with two kaons in 
the final state are therefore of particular experimental interest. Previous published upper 
limits on the branching fractions for the decays r _ — > £~Kg are of O(10~ 4 ) ||. There are 
no previous results for the decays r~ — > i^K^Kg. In this paper, we present the results of a 
search for the decays into one lepton and one or two Kg mesons, with the K s decaying into 
two charged pions. 

The data used in this analysis were collected using the CLEO detector ]7| from e + e~ 
collisions at the Cornell Electron Storage Ring (CESR) at a center-of-mass energy ~ 10.6 
GeV. The total integrated luminosity of the data sample is 13.9 fb^ 1 corresponding to the 
production of N TT = 1.27 x 10 7 r + r~ events. The CLEO detector is a general purpose 
spectrometer with excellent charged particle and shower energy detection. The momenta 
of charged particles are measured with three drift chambers between 5 and 90 cm from the 
e + e _ interaction point (IP), with a total of 67 layers. For ~63% of the data collected, the 
innermost tracking chamber was replaced by a three-layer silicon vertex detector M. The 
specific ionization (dE/dx) of charged particles is also measured in the main drift chamber. 
The tracking system is surrounded by a scintillation time-of-flight system and a CsI(Tl) 
calorimeter with 7800 crystals. These detector systems are installed inside a superconducting 
solenoidal magnet (1.5 T), surrounded by an iron return yoke instrumented with proportional 
tube chambers for muon identification. 

The r + r _ candidate events must contain four or six charged tracks with zero net charge. 
The polar angle 9 of each track with respect to the beam must satisfy | cos 9\ < 0.90. To reject 
beam-gas events, the distance of closest approach of each non-i^ track to the IP must be 
within 0.5 cm transverse to the beam and 5 cm along the beam direction. Photons are defined 
as energy clusters in the calorimeter with at least 60 MeV in the barrel (| cos^l < 0.80) or 
100 MeV in the endcap (0.80 < | cos 6 1 ! < 0.95). We further require every photon to be 
separated from the projection of any charged track by at least 30 cm unless its energy is 
greater than 300 MeV. In order to diminish QED background such as radiative Bhabha and 
/i-pair events with photon conversion, we require each event to have total energy less than 
95% of the center-of-mass energy. This requirement rejects most of the QED background, 
while incurring a small loss in detection efficiency. 

We divide each event into two hemispheres (signal and tag), one containing one charged 
track and the other containing three or five charged tracks, using the plane perpendicular to 
the thrust axis @. The thrust axis is calculated from both charged tracks and photons. The 
invariant mass of the tag hemisphere must be less than the r mass, M T = 1.777 GeV/c 2 [|TUJ. 
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The signal hemisphere must contain an electron or a muon and one or two K s mesons. The 
electron candidate must have shower energy to momentum ratio in the range, 0.85 < E/p < 
1.10, and when available, the specific ionization lost must be consistent with that expected 
for an electron. The muon candidate must penetrate at least three absorbtion lengths of iron. 
The Kg candidate is reconstructed in the 7r + 7r~ final state with a detached vertex, and the 
invariant mass must be within approximately three standard deviations of the nominal mass, 
485 < m wn < 510 MeV/c 2 , as determined from a signal Monte Carlo simulation (see below). 
In order to diminish radiative Bhabha and /x-pair events further, neither pion should be 
consistent with identification as an electron. 

Since there is no neutrino in the signal hemisphere while there is at least one neutrino 
undetected in the tag hemisphere, the missing momentum of the event must point toward 
the tag hemisphere, < cos 6 hag-missing < 1-0. In order to suppress the background from 
radiative Bhabha and /z-pair events, the direction of the missing momentum of the event is 
required to satisfy | cos missin g \ < 0.90. For the decay t~ — > e~K s , cos 6 ta g-missing is further 
restricted to be less than 0.99 to reduce the radiative backgrounds. This corresponds to the 
minimum ratio of background to detection efficiency The background is estimated from the 
sidebands in the invariant mass vs. total energy distribution of the decay candidates. 

To search for decay candidates, we select r candidates with invariant mass and total 
energy consistent with the expectations. The following kinematic variables are used to select 
the candidate events: 

AE = E — E beam 
AM = M — M T , 

where Ef, eam is the beam energy, and E and M are the reconstructed r energy and mass. 
The AE vs. AM distributions of the decay candidates in the data and signal Monte Carlo 
samples are shown in Figs. [I] and The center of the signal region is slightly shifted from 
zero in the AE vs. AM plane to account for initial state radiation and shower leakage. The 
signal region is defined as the area within three standard deviations (a) of the expectation 
for both kinematic variables, as determined from a signal Monte Carlo simulation. 

In the Monte Carlo simulation, one r lepton decays according to two- or three-body 
phase space for the mode of interest, and the other r lepton decays generically according 



to the KORALB /TAUOLA r event generator 11]. The phase space model is appropriate 
for an unpolarized tau. If the Lorentz structure of the neutrinoless decay is V-A, as in 
the Standard Model, correlations between the spins of the two r's in the event will lead to 
slightly higher detection efficiency than phase space, while V+A decays will lead to lower 



detection efficiency. The detector response is simulated using the GEANT program [[12 
The estimated detection efficiencies (e) (T^j are summarized in Table |. 

The upper limit on the branching fraction is related to the upper limit A on the number 
of signal events by 

B A 



S 



where B\ = (84.71 ± 0.13)% is the inclusive 1-prong branching fraction [l(J, B K o^ n+7T - = 
(68.61 ± 0.28)% is the branching fraction for K s to decay to two charged pions, and n 
is the number of Kg mesons in the final state. No candidate decays are observed, so we 
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FIG. 1: AE vs. AM distribution of the (a) data and (b) signal Monte Carlo sample for the decay 
t~ — > e _ -fC^; (c) and (d) show the corresponding distributions for r~ — > fi~Kg. The normalization 
of the signal Monte Carlo sample is arbitrary. The ellipses indicate the signal region (see text). 
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FIG. 2: AE vs. AM distribution of the (a) data and (b) signal Monte Carlo sample for the 
decay r~ — > e~KgKg; (c) and (d) show the corresponding distributions for r~ — > \x~ K^K^. The 
normalization of the signal Monte Carlo sample is arbitrary. The ellipses indicate the signal region 
(see text). 
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Mode e (%) ff(KT 7 ) (stat.) g(10~ 7 ) 

e~K° s 19.4 ±0.4 8.5 9.1 

(TK% 19.0 ±0.4 8.7 9.5 

e~K%K% 12.1 ±0.1 20 22 

^K° S K° S 8.0 ±0.1 30 34 

TABLE I: Summary of detection efficiency (with statistical uncertainty), 90% C.L. upper limits on 
the branching fraction with and without including systematic uncertainty. 



take A as 2.44 events in each mode at 90% confidence level, according to the frequentist 
method [f[3J| . The upper limits on the branching fractions with and without systematic 
uncertainties are shown in Table |. The systematic uncertainties include the t + t~ cross 
section (1%), luminosity (1%), track reconstruction efficiency (1% per charged track), K s 
detection efficiency (2% per Kg), lepton identification (1.5% for electron and 4% for muon), 
and the statistical uncertainties in the detection efficiencies due to limited Monte Carlo 
samples (1-2%). 

Black et al. [|15j have analyzed the constraints on the new physics scale for dimension- 
six fermionic effective operators involving r — fi mixing, motivated by the observed — v T 
oscillation. The most stringent lower limits from exotic heavy quarks and r decays on 
the physics scale of the operators involving quarks are ~10 TeV. The new upper limit on 
B(r~ — > fi~K s ) presented in this paper yields a lower limit of 17.3 and 18.2 TeV for the 
axial vector and pseudoscalar operators, respectively. 

In conclusion, we have searched for r decays involving Kg mesons that violate lepton 
flavor, but find no evidence for a signal. This results in improved upper limits for the decays 
t~ — > i~K s and first upper-limits for the decays r~ — > £~K S K S . The upper limits for the 



t — > £ K S K S final states are more stringent than those found previously for £ K + K [16 
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